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Sound driven gas bubbles in water can emit light pulses. This 
phenomenon is called sonoluminescence (SL). Two different phases of 
single bubble SL are known: diffusively stable and diffusively unstable 
SL. We present phase diagrams in the gas concentration vs forcing pres- 
sure state space and also in the ambient radius vs gas concentration and 
vs forcing pressure state spaces. These phase diagrams are based on the 
thresholds for energy focusing in the bubble and two kinds of instabili- 
ties, namely (i) shape instabilities and (ii) diffusive instabilities. Stable 
SL only occurs in a tiny parameter window of large forcing pressure 
amplitude Pa ~ 1.2 — 1.5atm and low gas concentration of less than 
0.4% of the saturation. The upper concentration threshold becomes 
smaller with increasing forcing. Our results quantitatively agree with 
experimental results of Putterman's UCLA group on argon, but not on 
air. However, air bubbles and other gas mixtures can also successfully 
be treated in this approach if in addition (iii) chemical instabilities are 
considered. - It is only the Rayleigh-Plesset ODE approximation of the 
bubble dynamics - extended in an adiabatic approximation to include 
mass diffusion effects - which allows us to explore considerable portions 
of parameter space. Therefore, we checked the adiabatic approximation 
by comparison with full numerical solution of the advection diffusion 
PDE and find good agreement. 
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I. INTRODUCTION 



A. The phenomenon 

A gas bubble levitated in a strong acoustic field 

P{t) = Pa cos tut (1.1) 

can emit bursts of light so intense as to be observable by the naked eye Here, 
Pa is the forcing pressure amplitude and uj/27r the frequency of the forcing field. This 
phenomenon is called single bubble sonoluminescence (SL). The light pulse is shorter 
than 50ps 0,0. Precise experiments by Putterman's group at UCLA have 
revealed many surprising and intriguing properties of sonoluminescing bubbles. SL 
only occurs in a narrow parameter range. The adjustable experimental parameters 
we focus on here are the forcing pressure amplitude Pa and the gas concentration 
Coo far from the bubble. Single bubble SL is found only for large (compared to the 
ambient pressure Pq) forcing pressure Pa ~ 1.2 — l.Satm and small (compared to the 
saturation cq) gas concentration Coo- 

We report the results on argon bubbles first. Two distinct phases of single bubble 
SL are known p: (i) Unstable SL occurs in the concentration range Cqo/cq = 6% — 
26%. This phase is characterized by an increase of the relative phase of light emission 
with respect to the driving pressure on the slow diffusive time scale ~ Is, followed by 
a rapid breakdown and another subsequent increase. The light intensity itself behaves 
in the same way and the bubble is reported to be dancing or jiggling |[l|J^. This state 
of SL is also unstable in the sense that often all of a sudden the bubble dies, (ii) 
Stable SL occurs in argon bubbles at very low gas concentrations Coo/cq ~ 0.4%. The 
diffusively stable state is characterized by the constancy of the relative phase of light 
emission over billions of cycles (days). The same is true for the light intensity. 

For air bubbles the same two phases exist, however, for gas concentrations about 
two orders of magnitude larger than for argon bubbles: Stable SL is observed for 
c'^/cq ~ 10 — 20% 0, unstable SL for even higher concentrations. 

Another controllable parameter is the temperature of the liquid. Upon decreasing 
the temperature of the ambient water from room temperature to slightly above freez- 
ing, the light intensity increases by two orders of magnitude ^ . Abrupt transitions in 
the light intensity with the liquid temperature are found for SL in nonaqueous fluids 
pi]| . Using different fluids (but the same gas species) also results in great differences 



in SL intensity ||TT|JTn]. 



B. A hydrodynamic approach 

The goal of this paper is to figure out which features of SL can be accounted for 
in a purely hydrodynamic approach ^ and to which extent they may reflect other, 
non-hydrodynamic effects, e.g. chemistry. 
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Our main result is the phase diagram figure |^ in the Pa — c^o phase space. It 
is obtained from hydrodynamic calculations of the bubble dynamics and the fiuid 
dynamical and diffusive processes outside the bubble. For given forcing pressure 
amplitude Pa and gas concentration Coo we predict with this diagram whether the 
bubble is in the diffusively unstable SL state, the diffusively stable SL state, or in no 
SL state at all. 




FIG. 1. Phase diagram in the Cqo/cq vs Pa/Po parameter space. The three phases 
represent stable SL, unstable SL, and no SL. For lucidity we do not draw the upper and the 
right borderline of the unstable SL regime towards a no SL regime as they are less precisely 
defined. 

Besides the forcing pressure amplitude Pa and the gas concentration Coo, two 
further parameters have to be addressed here: The frequency 00/2-71 of the forcing 
field and the ambient radius Rq of the bubble, i.e., the bubble radius at ambient 
normal conditions of Pq = latm and 293K. The frequency is set to a fixed value so 
that the forcing acoustic field ( |1.1D corresponds to a resonance of the container, in 
order that the bubble remains trapped in a pressure antinode. All of the analysis 
corresponds to = 27r ■ 2Q.5kHz as applied in Barber et al.'s experiment the 
period of the forcing field is thus T = 2Tr/uj = 38fis. 

The ambient radius Rq of the bubble is not an adjustable parameter but the 
system chooses Rq dynamically, as emphasized already by Putterman. The approach 
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followed here includes this (diffusive) dynamics, so that we can predict the ultimate 
ambient radius. Our results are consistent with Mie scattering radius measurements 
i- 

What are the necessary requirements for SL to occur? First, the bubble has 

T7[| . We have identified three types of 



to be stable towards shape oscillations 
shape instability: Two instabilities of a parametric type acting on relatively slow 
time scales of ~ 1 — lOO/xs. As discussed below, these instabilities are quite gentle 
and the bubble can survive them by pinching off micro bubbles. This pinch off causes 
the aforementioned break in the relative phase of light emission in the unstable SL 
state 0. The third type of shape instability however acts on a very short time scale 
< lO^^s. It is so violent that it ejects the bubble from the trapping sound field. We 
call it the Rayleigh- Taylor instability as it occurs when gas from inside the bubble is 
accelerated towards the fluid. 

What is the energy focusing process in the bubble? Many speculations abound 

^ (for multi bubble SL) and later Greenspan and 
suggest that shocks detach during the com- 



18,19 . Jarman 



in literature 

Nadim |2T]] and Wu and Roberts 
pression of the bubble and focus to the center of the bubble, thereby compressing 
the gas so strongly that light can be emitted, either by ionization and subsequent 
bremsstrahlung |2^J§] or by black body radiation [^. In ref. I^S] we suggest an 
alternate energy focusing mechanism. The idea is that the bubble acts as a driven 
acoustic resonator which switches on when the damping losses through viscosity and 
acoustic radiation are smaller than the energy input during the collapse. The acoustic 
energy accumulates and finally results in light pulses. The above mentioned abrupt 
transition in the SL intensity with increasing temperature |]10| can be accounted for 
within our theory, as well as the dependence of the light intensity on different liquids 
and gases. We refer to ref. for a detailed discussion. 

In this paper we take as criterion for energy focusing and the resulting light 
production that the Mach number of the bubble wall has to be larger than one 
. This criterion corresponds to the onset of SL within conventional shock theories 
20| , |22|pl| and within our alternate energy focusing mechanism Together with 



the requirement of bubble stability, it gives the boundary of the stable SL regime 
of Fig. ^ For argon bubbles the narrow parameter range where stable SL exists is 
in agreement with experiments of Barber et al. 0. A qualitative argument for this 
agreement was previously given by Lofstedt et al. 0. 



C. Chemical instabilities 

For air bubbles there are severe deviations between the hydro dynamically calcu- 
lated phase space diagram and the experimental measurements 0. The parameter 
regime where stable air bubbles should exist is close to that of argon; however, ex- 
perimentally stable SL is found for gas concentrations as large as /cq ~ 20% [Q. 
Because of this discrepancy, Lofstedt et al. hypothesize [0 a "yet unidentified mass 
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ejection mechanism" . In ref. |]25| we have suggested that this mechanism is chemicaL 
Indeed, when considering besides (i) shape instabilities and (ii) diffusive instabilities 
also (iii) chemical instabilities, our results can be extended to gas mixtures and are 
then in quantitative agreement with the UCLA experiments, as shown in detail in 
ref. The idea is that because of the high temperatures achieved in the bubble 
nitrogen and/or oxygen is destroyed and reacts to NO3 , NO^, and/or NH^ and only 
pure argon remains in the bubble. Thus for air which contains about 1% argon the 
gas concentrations in water have to be about two orders of magnitude higher than 
for pure argon. The central parameter is thus the argon (or inert gas) concentration 
^oo = qc^^^^^'^ in the dissolved gas. Here, q is the percentage of argon in the mixture; 
for air q = 0.01 = 1%. The nitrogen dissociation theory suggests that when adjusting 
q properly, no degassing is necessary any more p5|j26|] . 

In this paper we work out the basics of our hydrodynamic approach and restrict 
ourselves to pure argon bubbles for which no chemical instabilities (i.e., reactions) 
can occur. However, by considering the chemical instabilities properly [p5| , p6| , our 
results here can directly be extended to any gas mixture and are found to agree with 
the UCLA experiments. 



D. Necessary approximations 

How can we examine the huge multi-dimensional parameter space {Pa, Coo, and 
i?o)? Given that the dynamics involves time scales spanning eleven orders of mag- 
nitude (from the time scale of the light flash (< 50ps) to the diffusive time scale 
(~ Is)), it is necessary to make approximations in modeling the hydrodynamics of 
sonoluminescence. The full hydrodynamic problem involves solving the 3D Navier- 
Stokes equation both inside and outside the bubble, coupled with equations of heat 
transfer and gas transfer, accompanied by the correct boundary conditions at the in- 
terface and at infinity and the equations of state. Moreover, at least in principle the 
radiation fields need to be coupled to the fluid. This set of equations must be studied 
not only as a function of parameters but also over millions of oscillation periods of 
the bubble; it should be emphasized that the relevant question for sonoluminescence 
experiments is not the transient that occurs for the first few cycles but rather the 
nature of the long time limit. This complete formulation is both computationally and 
theoretically intractable. In order to make progress, the problem must be simplified. 
To date, two different avenues have been pursued: 



The first approximation was proposed long ago by Lord Rayleigh ||2^ and elabo- 



rated upon by Plesset |2^, Taylor Lauterborn Prosperetti pT| , p!6| and others 



|3^ , in the context of studies of cavitation. The idea is to consider the bubble as 



a perfectly spherical cavity, with the pressure inside the bubble having no spatial 
variations. The temporal variation of the pressure follows from an equation of state. 
In this approach the full dynamics is reduced to the Rayleigh-Plesset ODE [03. This 



formulation allows very long time calculations of the bubble dynamics, but it com- 
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pletely ignores the dynamics inside the bubble producing the hght 
Strube 
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Later on, Plesset 

and Prosperetti |jl6| extended this type of approach to deal with 



shape oscillations while Epstein and Plesset , Eller and Crum fS^JT^ , Crum and 
Cordry and finally Fyrillas and Szeri and Lofstedt et al. included diffusive 
effects. We call this approach the RP-SL-bubble approach. Clearly, shock formation 

[23[] inside the bubble will modify 



p0| , |21|| or the building up of the acoustic waves 
the dynamics of R{t). But it is our belief that the results of this paper are robust 
towards these changes and it is only within this RP-SL bubble approach that the 
exploration of the full SL parameter space and the calculation of phase diagrams are 
currently manageable. Full numerical simulations as in [^,0 are by far numerically 
too expensive to do such an analysis. 

The second type of approximation traditionally made is complementary to the 
first, and focuses on the interior of the bubble and above mentioned shock formation 
processes |^J2^,|5S|,^ . The spherically symmetric gas dynamics equations are solved 
inside the bubble and coupled to the Rayleigh-Plesset equation. Simplifications are 
typically employed in modeling the gas dynamics, for example neglecting heat and 
viscous dissipation. These calculations can only be carried out for a few oscillation 
periods, and thus are not able to resolve cumulative effects building up over many 
oscillations. 



E. The Rayleigh-Plesset equation 

The Rayleigh-Plesset (RP) equation |^,27,30|, on which the entire analysis of this 
paper is based, describes the dynamics of the bubble radius, 

RR +^R' = — {p{R, t) - P{t) - Po) 

pwCiu at n pyjJri 

Typical parameters for an argon bubble in water are the surface tension a = 
0.073kg/ s"^, the water viscosity u = 10~^m^/s, density = 1000kg / , and speed of 
sound = 1481m/s. The driving frequency of the acoustic field is (^/2'7r = 26.5kHz 
and the external pressure Pq = latm. We assume that the pressure inside the bubble 
varies according to 

where h = i?o/8-86 is the hard core van der Waals radius. The exponent 7 is the 



effective polytropic exponent of the gas. Plesset and Prosperetti [|3T| calculated how 
it depends on the (thermal) Peclet number Pe = Rqu/k which gives the ratio be- 
tween the bubble length scale Rq (which we take as ~ 5/im for the estimates in this 
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paragraph) and the thermal diffusion length ^JhJuj. The thermal diffusivity k for 
argon is k ^ 2 ■ 10~^m^/s, which yields Pe ~ 0.2 and according to fig. 1 of |HT|], the 
effective polytropic exponent 7 = 1. As discussed below, the RP equation contains 
much smaller time scales than lo^^ . One could therefore argue that these smaller 
time scales may enter into the calculation of Pe, so that the frequency lo should be 
replaced by |P|/P. This estimate would lead to Pe(t) as large as 10^ at instants 
of rapid bubble wall movement which implies 7 ~ 5/3 for argon. However, since 
Pe(t) ^ 1 only holds in very small time intervals ~ Iras, the global dynamics are not 
affected by setting the effective polytropic exponent 7 = 1 uniformly in time. Note 
that with 7 = 1 eq. ( p..3| ) should not be thought of as equation of state but rather as a 
process equation parametrizing the isothermal conditions at the bubble wall, induced 
by the large heat capacity of water. The choice of 7 = 1 is confirmed by the full 



numerical simulations of Vuong and Szeri PS| and by the approximation of Kamath 



et al. Note that, as a consequence, there are heat fluxes back and forth across 
the bubble wall. 

The radius P(t) corresponding to the forcing pressure ( |1.1|) with P^ = l.lSatm 
is shown in figure j^. Four time scales are hidden in the RP equation: The pe- 
riod T = 38/is of the external forcing P(t) (figure ^), the intrinsic frequency 

^37Po/(p^Po)/27r ~ (1.8;us)~^ of the oscillating bubble which is the frequency of 
the afterbounces, the time scale of viscous damping Pq/'^ ~ 25/is, and the duration 
~ 0.1 — Ins of the bubble collapse, estimated in H. A fifth time scale, determined 



by the surface tension yPop^/cx ~ 1.3/iS is only important for large a or small P; for 
typical Po, it only slightly changes the intrinsic time scale which is of the same order 
of magnitude. 

Also diffusive processes can be understood within the RP approach to the SL 
bubble [0,^,^]. At first sight this is surprising because there is no diffusive time 



scale in the RP equation. For a qualitative understanding why this works nevertheless 
look at the bubble radius R{t) (figure |^), resulting from ( |1.2| ). For large R(t) the 
pressure inside the bubble will be low and gas diffuses into the bubble (rectified 
diffusion). For small P(t), on the other hand, the bubble will shrink because of 



the enormous pressure inside |43|. This concept was made quantitative by Fyrillas 



and Szeri and Lofstedt et al. 0. The main idea is a separation of the slow 



diffusive time scale from all time scales in ( |1.2| ). We therefore call this approach the 
adiabatic approximation of diffusion. As we will see the balance between growth and 
shrinking is very delicate. In fig. |^c we present the ambient radius (calculated in 
section 0) Po(t) which represents the mass m = AnRlpo/S of the bubble, with the 
ambient density po ~ l.Qkg/m^ for argon. It corresponds to the R{t) curve in fig. 
||b. The exchange processes between the bubble and the liquid can be very violent. 
In fig. ^ the bubble's mass increases to 100.5% of its initial value and decreases to 
99.8% thereof while after one full cycle it again takes its initial value. These are 
representative values for argon bubbles near the onset of the SL regime; for larger 
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forcing, the exchange processes become even more violent. 
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FIG. 2. (a) Forcing pressure P{t) = PaCOSLot, Pa = l.lSatm for two cycles and the 
corresponding (b) R{t) and (c) Ro{t). The bubble is near an equilibrium state. The gas 
concentration is Cqo/cq = 0.035. 



F. Organization of the paper 

The paper is organized as follows. In section II we analyze the bubble stabil- 
ity with respect to shape oscillations. We then give as necessary criterion for light 
emission energy focusing in the bubble, either through a shock or through acous- 
tic resonance (section |IT|). Section |I^ constitutes the main part of the paper. We 
calculate the diffusive instability with Fyrillas and Szeri's and Lofstedt et al.'s 
@] adiabatic approximation. The main result are phase diagrams in the Coo — Pa, 
the Ro — Pa, and the Rq — Coo parameter spaces. Only in a very small parameter 
domain does stable SL occur. To analyze unstable SL we calculate growth rates and 
compare them to Putterman's measurements of diffusively unstable SL argon bubbles 
1^. Concentration profiles from a full numerical solution of the advection diffusion 
PDE are presented in section 0, where we also check the validity of the adiabatic 
approximation by comparison with full numerical solutions of diffusive dynamics. We 
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find experimentally undetectable discrepancies which vanish in the Schmidt number 
Sc oo limit. Section presents conclusions. 



II. SHAPE STABILITY 



For sonoluminescence to occur and for the bubble to remain oscillating for bil- 
lions of cycles, the bubble must be stable to shape oscillations. First, following the 
pioneering work of Plesset [|l^], Strube |T^, and Prosperetti [0, we derive equations 



Strube fll 

for the deviations of the bubble from a spherical shape, and then proceed to analyze 
them. 



A. Dynamical equations 



We focus on the stability of the radial solution R{t). Consider a small distortion 
of the spherical interface R{t), 

R{t) + an{t)Yn{0,<P), 

where Yn is a spherical harmonic of degree n. The goal is to determine the dynamics 
a„(t) for each mode. Plesset's [|12| derivation follows the same spirit as the derivation 



of the Rayleigh-Plesset equation. The potential flow outside the bubble is constructed 
to satisfy the boundary condition that the velocity at the bubble wall is i? + hnYn- 
This potential is then used in Bernoulli's law to determine the pressure in the liquid at 
the bubble wall. Applying the pressure jump condition across the interface yields the 
Rayleigh-Plesset equation for R{t) as well as a dynamical equation for the distortion 
amplitude a„(t). 



3i?. 

ri 



[n 



1 



R PnCr 



R Pv^R^ 



0, 



(2.1) 



where /?„ = (n — l)(n + l){n + 2). However, viscous effects have been neglected in 
Plesset's derivation. 

The intrinsic difficulty 
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Viscosity was later taken into account by Prosperetti 
in its consideration is that viscous stresses produce vorticity in the neighborhood of 
the bubble wall [^. In principle, vorticity spreads both by convective and by diffusive 
processes all over the fluid and the problem becomes nonlocal. However, for small 
viscosity the generated vorticity will be more or less localized and we can introduce 
a bubble boundary layer approximation of the nonlocal equations which we do in the 
next subsection. 

Here, we give the dynamics of the nonlocal problem, closely following Prosperetti 
p!6| . It is advantageous to decompose the vorticity field in the fiuid in a poloidal and 
a toroidal part, which are conveniently represented by scalar fields S{r, t) and T(r, t), 
respectively. 
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a; = V X V X [5(r,t)F-(^,</))er] + V x [T(r,t)y-(^, 



(2.2) 



Only the latter, T{r,t), contributes to the long term dynamics of the bubble. Its 
dynamics is given by the PDE 



d,T{r, t) + R'Rdr i \t) = udlT - ^^I^il^T, 



by the nonlocal boundary condition at R{t), 

POO 2 

T{R, t) + 2i?"-i / s-^'Tis, t)ds = 

JR n + 1 



[n + 2)dn-in-l)^R 



(2.3) 



(2.4) 



and by the boundary condition at infinity, T{oo, t) = 0. At t = the fiuid is assumed 
to be at rest. The collapse R{t) of the bubble transports vorticity u: into the fiuid, 
see eq. ( |2.3| ). 

Once created, the vorticity acts back on the dynamics of an{t). These indirect 
viscous corrections together with the direct ones modify eq. ( |2.1| ) to yield 



with 



dn + Bn{t)an - v4„(t)a„ + n{n + l)(n + 2)—T{R, t) 

R'^ 



s 



-j T(s,t)c/s = 



(2.5) 



AJt) = {n- 1)- - + 2/3„z/4- 

^ ^ ^ ' R p^R^ ^ R^ 



(2.6) 



Bn{t) 



3R 
~R 



V 

R^' 



(2.7) 



B. Boundary layer approximation 



For an exact stability analysis the coupled eqs. ( |2.3| ) - ( |2.5| ) together with the RP 
equation (|l]2|) must be solved. However, considerable vorticity is only to be expected 
in a small boundary layer of thickness 6 around the bubble. 

Within this boundary layer approximation the space integrals in ( p. 41) and 
can be approximated by the (integrand at R)x6. The integral in ( |2.4| ) thus is 
R~^T{R,t)6 while the one in (|2.5|) vanishes. We obtain 

dn + Bn{t)an - An{t)an = (2. 



with 
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R (3n(r 2uR 



1 + 2 



(2.9) 



3R 2u 
~R^^ 



n{n + 2f 
1 + 2 ' 



(2.10) 



The viscous contribution to A„(t) is not important and only causes a tiny shift, 
as the ratio between the third and the second term of the rhs in ( p.9| ) is typically 
ppwRQUj/a ~ 10^^. However, in ( p.lO|) it introduces a damping rate 



n(n + 2)2 
1 + 2| J 



(2.11) 



acting on the shape oscillations. That only the second term in (|2.10 ) contributes to 



the damping rate of the oscillator can formally best be seen after the substitution ||T2 



bn{t) oc (i?(t))^/^a„(t). Physically this is not surprising, as the first term is a mere 
consequence of spherical geometry. Two physical effects contribute to the damping 
rate: (i) Stabilizing, local damping by viscous dissipation. In this limit the viscous 
boundary layer around the bubble vanishes {6 = 0) and the damping rate becomes 
^nit) = 2i'{n + 2)(2n + l)/i?^ > 0. (ii) The movement of eddies around the bubble, 
generated by the shape oscillations itself. With increasing boundary layer thickness 
6 (i.e., with increasing viscosity) this destabilizing effect becomes stronger. 

How to approximate the thickness 6 of the boundary layer when we have non- 
vanishing vorticity? For large bubbles ^ 5 it is set by the diffusive length scale 



^Jv/uj = 2.5/zm in equation ( p.3| ) p5 



As typical frequency scale we choose the forcing 
frequency u. Higher frequencies are of course also present in the RP dynamics, but 
a Fourier analysis of the R(t) signal shows that the forcing frequency is dominant. 
Because of the angular contribution to the dissipation (the second term on the rhs 
of eq. (|2.3|) ) we also expect a slight dependence on the spherical mode n which we 
neglect here. 

For small bubbles -R <^ 5 we do not expect the boundary layer around the bubble 



to be larger than the bubble itself. We thus have to introduce a cutoff EH] . We choose 



6 



mm 



V R 

uj ' 2n 



(2.12) 



The n-dependence of the cutoff can be understood from the quasi static limit which 
holds for small bubbles, as for small bubbles the bubble dynamics is strongly damped 
by viscosity and surface tension and R{t) does not change much. In this quasi static 
limit the Ihs of ( [2.3[ ) vanishes and T(r) = T{R){r / R)~^ is the static solution. It 
decays to half its boundary value T{R) at r = 2^/"/?. Thus 6 = i?(2i/" - 1) ^ 
R\n2/n ^ R/{2n) as in (^121) . More precisely, for T(r) = T{R){r/R)-'' we can 
calculate T{R) from eq. ( p.4| ), and obtain essentially the same S. The exact values 
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of our results depend on details of the cutoff ( |2.12| ). However, the general features of 
the solution are invariant. 

With the approximation ( p.l2| ) we can Taylor-expand ( p.9| ) and ( |2.10| ) and finally 
obtain as approximate dynamical equation for Qnit) equation ( p.8|) with 



R I3n(y 2uR 



(n - l)(n + 2) + 2n(n + 2){n- 1) 



R 



(2.13) 



Bn{t) 



3R 2u 

'r^r^ 



(n + 2)(2n+ 1) - 2n(n + 2)^ 



R 



(2.14) 



Our results are based on these equations. 



C. Ray leigh- Taylor, afterbounce, and parametric instabilities 

Three types of shape instabilities are hidden inside these equations. We call them 
the Rayleigh- Taylor instability, the afterbounce instability, and parametric instability, 
for reasons which will become clear later. They are distinguished by the widely 
different time scales over which they act. The transition between these instabilities 
is often gradual rather than abrupt. Nevertheless, we think that our classification is 
physically important as the difference in their time scales results in a difference of the 
typical velocities of the bubble fragments after the shape instability has destroyed 
the bubble. We estimate this velocity as typical length scale ~ 1/xm of a collapsed 
bubble divided by the typical time scale of the pinch off. 

The goal of this subsection is to find criteria for the occurrence of the three shape 
instabilities. Their nature becomes clear from figure ^, where we show the dynamics 
of a2{t) (which is the most unstable spherical mode in the parameter range discussed 
in this paper, so we restrict ourselves to it) for the three different regimes of instabil- 
ity. We normalize 02 (t) to the current bubble radius R{t) to get a measure of bubble 
distortion. The upper part displays the dynamics of the distortion amplitude 02 (t), 
normalized to R{t), in a Rayleigh- Taylor unstable, parametrically stable case, ob- 



tained from a numerical solution of eqs. ( |1.2|) , ( p.8|) , (|2.13|) , and (|2.14|) , with addition 
of small amplitude noise. The middle part shows the dynamics of the distortion in 
the regime of the afterbounce instability, and the lower part typifies the dynamics 
of a parametrically unstable, Rayleigh- Taylor stable bubble. Clearly, very different 
time scales are responsible for the distortion of the spherical shape of the bubble. Of 
course there are regimes in the Rq — Pa parameter space where the bubble is stable 
(or unstable) towards two or even all three instabilities. 

First we focus on the Rayleigh- Taylor instability, occurring near the minimum 
bubble radius when the gas accelerates into the fluid. The strongest destabilization 
occurs just after the bubble radius reaches its minimum. The acceleration of gas 
towards the fluid during this time is enormous, motivating the name Rayleigh- Taylor 
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instability. Closer analysis of this section of bubble movement reveals that the 
time scale of the Rayleigh- Taylor instability can be estimated by the expansion time 
scale of R(t) just after the collapse which is tRT ~ h/cw with the van der Waals hard 
core radius h = i?o/8.86. Thus, t^T ~ 10~^ — 10~^°s, which is confirmed by the 
numerical results. In order to take into account microscopic fluctuations we added a 
random displacement of size ~ O.lnm to the distortion 02 (^) after each integration 
time step. 

For lower Pa the destabilization during the violent bubble collapse may not be 
strong enough to immediately overwhelm the bubble. But as seen from figure ^ 
further periods of destabilization occur during the afterbounces. As pointed out 
above, in the afterbounce regime R{t) oscillates on the bubble's intrinsic time scale 



tj ~ \J PujRq/Pq ~ This is too fast for viscous effects to smooth out the shape 
distortions, so after a couple of afterbounces the bubble may be overwhelmed. This 
type of instability shows features of parametric instability, however, the afterbounces 
are not strictly periodic. As an approximate criterion for this afterbounce instability 
we give that microscopic fluctuations can overwhelm the bubble within one period T, 

The transition between Rayleigh- Taylor and afterbounce instabilities is illustrated 
in figure The R time series shows violent behavior at the main bubble collapse 
and the first afterbounces. Then, the behavior of a2{t)/R{t) becomes oscillatory and 
locks into the periodicity of the R{t) afterbounces with twice their period (the same is 
true for 02 (t) itself). This is to be expected from a Mathieu type instability, which is 
most effective for a driving with twice the intrinsic frequency of the driven equation. 

Figure ^ depicts a phase diagram of shape instabilities as a function of the ambient 
bubble radius Rq and the forcing pressure amplitude Pa for n = 2. The dashed line 
gives the combined stability threshold for Rayleigh- Taylor and afterbounce instabili- 
ties (the "fast" instabilities with time scales <^ T). The global features of the phase 
diagram are easily understood. Small bubbles are more stable than large ones thanks 
to surface tension and viscosity, as the second term in ( p.l4p becomes dominant for 



small radii R. Evidently, weakly forced bubbles are more stable than strongly forced 
ones. 

A pure parametric shape instability acts on the much longer time scale of the 
forcing T ^ 38fis. It corresponds to a net growth of a nonspherical perturbation 
over one oscillation period, so that after many periods perturbations overwhelm the 
bubble. The time scale of the parametric instability tpj is thus many forcing periods 
T = 211/00. However, as the finally resulting pinchoff occurs during afterbounces, its 
time scale is the same as for the afterbounce instability, i.e., the intrinsic time scale 
of the bubble motion tj ~ 

In the relevant parameter regime for the parametric instability R{t) and thus also 
An{t) and -B„(t) are strictly periodic in time with frequency 1/T. Thus Eq. ( p7^ ) is 
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an ODE of Hill's type and the parametric instability can be rigorously analyzed. It 
occurs whenever the magnitude of the maximal eigenvalue of the Floquet transition 
matrix Fn{T) of eq. ( |2.8| ) is larger than one. The Floquet transition matrix Fn{T) is 
defined by 



Fn{T) " . (2.16) 





By numerically computing the eigenvalues of the Floquet transition matrix we mapped 
out the phase diagram of stability. Fig. |^ shows the stable and unstable domains in 
the Rq — Pa parameter space. 

In the SL parameter range of ~ 1.2 to l.Satm the bubble becomes parametri- 
cally unstable at about Rq^ ^ 4 — 5fim. This number is not to be understood as a 
prediction of the exact value, as within our approximations we can only predict the 
order of magnitude and trends. 

For smaller Pa the threshold for instability Rq^ does depend on the forcing pres- 



sure. We discussed phase diagrams in those regimes in [jT^ and also showed that in 
the small forcing limit eq. ( ^.81 ) reduces to a Mathieu equation. 

That Rq^ does not significantly depend on Pa for large Pa can be understood 



from the dynamics of R{t) and a„(t) and from equation (p.l4|) . If R(t) is small, 
the second term in ( p.l4| ) dominates and stabilizes an{t). For small Pa the minimal 
radius Rmin = mint{R(t)) still decreases with increasing P^. But after the van der 
Waals hard core radius h = Rq/8.86 has once been reached for large enough Pa, Rmin 
becomes Pa independent . 

All calculation have been performed for the viscosity of water u = 10~^m^/s. Of 
course Rq^ and the other thresholds strongly depend on u. E.g., for a viscosity five 
times that of water we have Rq^ ~ 10/xm, but we won't discuss this dependence here. 



D. After the shape instabihties 

All types of shape instabilities result in the pinching off of micro bubbles. In 



other experimental configurations such as the Faraday experiment |2^, parametric 
instability can saturate at a finite amplitude. This is also possible for larger bubbles 
which are driven with small forcing pressure Pa < latm • However, we believe 
the nonlinear saturation is unlikely in the present experiment because the bubble size 
changes by two orders of magnitude during a single oscillation period. Saturation 
would require that the amplitude a„ is much smaller than the minimum radius Rmin] 
however, since the bubble spends most of the cycle with R ^ Rmin, the nonlinearities 
mainly act at larger radii. 

The major question, therefore, is what happens after the pinch off of a micro 
bubble. Will the bubble remain trapped at the anti-node of the pressure field, or will 
it escape from the system? 
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The force holding the bubble in the trap is the so-called Bjerknes force given by 
11,11 



F = -^rfVP. (2.17) 

Both R{t) and VP(t) are oscillating with time. The combination of both will lead 
to an effective, period averaged force (F) which pushes bubbles smaller than the 
resonance radius (~ lOO/im here) to the pressure antinode |^. The bubbles are thus 
trapped. 

During the collapse the forcing pressure is positive and the force ( p. 171) repulsive. 
If micro bubbles pinch off at this instant, will they and the recoiled bubble be fast 
enough to escape from the node before the force becomes attractive again? If so, they 
must travel a quarter of the acoustic wavelength in time T/2, thus their velocity 

must be of order ~ ~ lO^m/s. Assuming a typical length scale of Rq ~ l/im for the 
micro bubble and the remainder, we obtain as critical time scale of the collapse 10~^s: 
If the collapse is faster, the bubble can not survive in the trap, if it is slower, it is 
likely to survive as a smaller bubble. From the estimates of the relevant collapse time 
scale in the last subsection we conclude that the bubble will survive the afterbounce 
and the parametric instability where the pinch off occurs on the intrinsic time scale 
tj ~ Ifis. The pinched off micro bubble from those shape instabilities will dissolve 
by diffusion (see below), but the remainder of the bubble can survive and grow by 
rectified diffusion. 

The bubble fragments can, however, be ejected from the trap after a RT instability 
{tRT ~ 10~^s). To get a more stringent criterion for ejection, we computed the 
maximal velocity \R + a2\- For a bubble split apart by shape instabilities, this will 
also give the typical fragment velocity. The outward velocity is dominated by the 02 
term; R does not exceed O.lCm ||4^. Therefore, we assume that the fragments will 
escape if 

max (\^]>1. (2.18) 

{t'\t<t'<t+T} \ / ~ 

Figure || depicts the Ihs of this inequality as a function of Pa and Rq. The per- 
turbation speed exceeds the sound speed in two clearly distinct regions of parameter 
space: at large radii, where the parametric instability leads to high velocity, and at 
small radii and high pressures, where Cy^ is reached during Rayleigh- Taylor instabil- 
ity. Of course, only the latter region will determine the escape of the bubble, because 
parametrically unstable bubbles shed micro bubbles long before these high velocities 
are reached. We extract a |a2| = c^, isoline from this graph and add it to the phase 
diagram in figure ^ (dotted line). In accord with the estimate for t^T in the pre- 
ceding section, bubbles with smaller R^ have faster instability time scales and thus 
higher values of |d2|. However, for very small Ro, the bubble dynamics is stabilized 
by surface tension, and the perturbation velocity drops again. The exact position of 
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this line is not meant to be quantitative, because it is extracted from dynamical data 
for bubble velocities where the RP equation is not a good approximation to bubble 
dynamics any more. 



III. ENERGY FOCUSING MECHANISM 



Although this paper is primarily concerned with stability constraints on a bub- 
ble obeying the Rayleigh Plesset equation, in order to relate the calculations to the 
sonoluminescence experiments it is necessary to adopt an onset criterion for the light 
emission. The light production is triggered by a hydrodynamic energy focusing mech- 
anism, which concentrates the input energy enough so that light is produced. The 
mechanism through which the focused energy produces light is outside the scope of 
this paper; many suggestions such as bremsstrahlung or black body radiation are 
present in literature p2|,|5|. 

Two theories of hydrodynamic energy focusing have been proposed: The original 
theory was that during the collapse of the bubble shocks detach from the gas-water 



interface and focus to the bubble's center pl]j2^ . The motion of a focusing shock 
is described by Guderley's similarity solution to the hydrodynamic equations pOj^ , 
which dictates that the temperature at the shock diverges as T cx Rj^, where Rg is 
the distance of the shock from the origin and j9 ~ 1 is an irrational scaling exponent. 
The amount of energy focusing in this theory is determined by the minimum distance 
Rs,min to which the shock approaches the origin. 



The second theory of energy focusing PB[ posits that the bubble is an "acoustic 



resonator", and that acoustic energy builds up in the bubble over many oscillation 
periods. Within this picture the amount of energy focusing is set by the total stored 
energy in the bubble. 

The crucial issue for the present paper is the onset criterion for how strongly the 
bubble must be forced for significant energy focusing to occur. It turns out that both 
shocks and collective energy buildup require the same onset criterion: the inward 
bubble wall velocity must be of order of the speed of sound in the gas, 

M = — >1. (3.1) 



Cgas 



Here, M is the Mach number, and 



2 p R 

'^9as = 7-^^Zrh^ (3.2) 

is the (squared) velocity of sound in a van der Waals gas. For the shock wave theory 
the criterion ( |3.1| ) is motivated by the fact that shocks can only detach from the 
bubble wall if its velocity is supersonic p^,|6|. For the acoustic resonator theory 



the criterion (|3.1|) holds as resonators can only be supplied with energy if they 
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are forced on a time scale comparable to the intrinsic time scale. Here, the forcing 
time scale is R/\R\ and the intrinsic acoustic resonator time scale is R/cgas- Thus 
\R\/cgas ~ 1- From above discussion of the time scales in the RP equation it follows 
that all time scales between R/\R\ and T are present in the R{t) dynamics; thus (|3.1| ) 
is the correct criterion for energy focusing. 

For large enough forcing the energy focusing criterion (|3.1|) is fulfilled once per 
cycle, namely shortly before the bubble achieves its minimum radius which in the 
relevant Pa domain is very close to the hard core radius h. This is where the light 
pulse is emitted |ll|,|^j36[] which gives support to the criterion (|3.1| ). In figure || we 
plotted the threshold for the M = 1 criterion in the Pa — Rq parameter domain. 

Let us check two further conditions which should be fulfilled within the shock 
wave theory. First, another requirement besides ( |3.1| ) is that the shock must persist. 
Dissipative mechanisms inside the bubble must therefore be weak. A measure for 
the relative strength of inertial and dissipative mechanisms is the Reynolds number 
Re. Viscous effects dominate for small Re. As crossover for nonlinear effects such as 
shocks to take over we take 112^ 



Re = ^ > 10. (3.3) 



^gas 



The kinematic viscosity of argon is Ugas = 11 ■ lO^^m^/s. We neglect its temperature 
dependence. In the center of the bubble this is a poor approximation and Ugas will be 
lower. The Re > 10 criterion will then be fulfilled earlier. However, ( |3.1|) is the more 
stringent criterion anyhow as seen from figure 

Second, we have to compare the thickness of a shock with the bubble's size. The 
thickness of a shock is of the same order of magnitude as the mean free path I of a gas 
molecule [^|. We have / ~ V/{Na()) where V is the bubble volume, N the number of 



particles in the bubble, and cxo ~ 10 ^^m^ the collision cross section of argon atoms. 
Thus 

With Iq ~ 10~^m we obtain / ~ lO^^^m during the strongest compression, i.e., a very 
sharply defined shock with a width <^ i? in spite of viscosity. 

We now come back to figure |^. That plot summarizes the criteria we suggest 
to be necessary for SL to occur: (i) Bubble wall Mach number M > 1 to ensure 
energy focusing to reach the high temperatures necessary for SL. (ii) Short time scale 
shape stability (Rayleigh-Tarylor and afterbounce) and (iii) parametric stability, (iv) 
Finally, the perturbation speed must not exceed c^, in order to keep the bubble or 
its fragments trapped in the sound field. 

There is only a small domain in parameter space where the bubble fulfills all four 
criteria. This domain is shaded in fig. It is this domain where we expect SL to 
be possible (within our RP-SL approach). However, up to now no statement on the 
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diffusive stability has been made. We will address this subject in the next sections 
and find that for low enough gas concentration the bubble in the shaded domain is 
also diffusively stable. 

IV. DIFFUSIVE STABILITY 

Two types of SL in argon bubbles have been observed 0: For large argon concen- 
trations Poo = 200mmHg or SOmmHg the SL bubble is diffusively unstable whereas 
for low concentrations Poo = 3mmHg the bubble is diffusively stable and the relative 
phase of light emission stays constant for days, see figure p!T| c for the experimental 
result. With the ambient pressure Pq = 760mmHg = latm these three gas con- 
centrations translate into relative concentrations of Cqo/cq = 0.26, 0.06, and 0.004, 
respectively, where cq = 0.061kg /m^ is the saturation (mass) concentration of argon 
in water for room temperature. 

In this section we set out to quantitatively understand the difference between 
the high and low concentration. Our goal is to calculate a phase diagram in the 
parameter space of the two experimental control parameters concentration Coo/cq and 
forcing Pq/Pq. For given concentration and given forcing we will thus be able to 
predict which of the three phases "diffusively stable SL", "diffusively unstable SL", 
and "no SL" will be realized. The third parameter, the ambient bubble radius Pq is 
not at the experimenter's disposal but the system will choose Rq itself. The ambient 
radius will follow from our analysis. First, we will present phase diagrams in the 
ambient radius - concentration and ambient radius - forcing pressure phase space. 

We again stress that it is the application of the adiabatic approximation of the 
diffusive problem p7|,|7[] which allows for the exploration of the whole 3D phase space 
-Ro ~ Coo — Pa- Strictly speaking, it only holds in the limit of zero diffusion constant 
D 0. However, in the next section we will show that for the physical diffusion 
constant D = Dat = 2 ■ 10~^m^/s the deviations between the exact solution and the 
adiabatic approximation are tiny, so that we can apply it here. 

A. Formulation of the diffusive problem 

Assuming spherical symmetry, the mass concentration of gas c(r, t) dissolved in 
the liquid at distance r > R(t) from the center of the bubble obeys the advection 
diffusion equation 

r'^dtc + R^RdrC = Ddrir^drc). (4.1) 
As boundary condition at the bubble wall we assume Henry's law 

c{R,t) = Cop{R,t)/Po. (4.2) 
The concentration at r — ^ oo is given by Coo, 
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c(00, t) = Co 



(4.3) 



The concentration gradient at the moving boundary gives the mass loss/gain of the 
bubble 



m = AuR^DdrC 



(4.4) 



The bubble is driven by the RP equation ( [1.21 ). Together with the initial conditions 
R(t = 0) = -Ro, R(t = 0) = and c{r,t = 0) = Coo this set of equations completely 
defines the problem. 



With the transformation to the Lagrangian coordinate ^5 

h{r,t) = ^{r'-R'{t)) 
the advection diffusion PDE ( |4.1| ) simplifies to the diffusive equation 

,4/3 



dtc{h,t)=Ddh 



{3h + R%t)) dhc{h,t) 



(4.5) 



(4.6) 



Eq. (|4.6|) can still not be solved analytically. How it can best be solved numerically 
is sketched in appendix 0. 



B. Adiabatic approximation 



The main idea of Fyrillas and Szeri |]3^ and Lofstedt et al. [0 is to treat the 



diffusive PDE by the method of separation of time scales |52]. They split the con- 
centration field in an oscillatory part Cosci'f', t) changing on the (fast) time scale T of 
the driving field and a smooth part Csmo{r,t) changing on a slow diffusive time scale 

c(r, t) = Coscir, t) + Csmoir, t). (4.7) 

This approach can be thought of as having introduced an "adiabatic" or slow time 
t. The smooth profile Csmo only depends on the adiabatic time, Csmo{r,i). In the 
PDE for Csmo{r,t) the fast time scale ~ T is averaged out. We define = Rq/D 
as diffusive time scale. Then the Schmidt number Sc = 2tttd/T is a measure of the 
quality of time scales separation. Sc ^ oo or D ^ means perfect separation. 
It turns out to be useful to introduce weighted time averages, 



j ^f{t)R\t)d t 



which may still depend on the adiabatic time t. Here, one only needs (.)jq and (.)j4. 
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The main result of ref. is that in the asymptotic hmit t 
profile Csmo{h,i) converges to 



,{h) = Coo + 



co- 



(Pit)) 



tA 



rh dh' 

JO ((3h'+/?3(t))4/3^^^^ 

roo dft/ 

•^0 ((3/l'+i?3W)4/3)^_^ 



oo the smooth 



(4.9) 



From 



the adiabatic growth of the bubble can be calculated as 





'c^ (P{t)>t,4©' 




Dco 


CO 





dh' 



{{3h'+RHt))i/3)^^ 



(4.10) 



The determination of the adiabatic growth rate has thus been reduced to solving the 
RP ODE ( |1.2| ) for R{t), calculating time averages (.)^ • of functions of R{t), and the 
solution of an space integral. We thus understand the adiabatic approximation as 
being in the spirit of the RP approach. 



C. Equilibrium points 



We now apply eq. ( [4.10|) to the acoustically driven argon bubble 0. Let us 



disregard the shape instabilities discussed in section || for the time being. 

The ambient bubble radius is in equilibrium (within the adiabatic approximation) 

if 



Coo (p(^))/4 



Co Po 



(4.11) 



Note that within the adiabatic approximation the condition does not depend on the 
diffusion coefficient D. However, the smaller D, the better the adiabatic approxima- 
tion holds. The equilibrium is stable, if 

(4.12) 

is positive. Scaling laws for (p)i4 and (3 are discussed in the next subsection and in 
ref. here we concentrate on the physical consequences. 

In figure |^ we plot {p{t))^^ as a function of Rq for various forcing pressure am- 
plitudes Pa- We first focus on small Pa ~ O.Satm and high gas concentration of 
Coo/cq ~ 0.7. There exists an unstable equilibrium at Rq ~ 6/im. Smaller bubbles 
shrink and finally dissolve, larger bubbles grow by rectified diffusion. For larger Pa 
the average {p{t))^^ (Rq) behaves quite differently in the small Rq regime. It starts to 
show characteristic wiggles, which can also be seen in Rmax{Ro)- Here, Rmax is the 
maximal radius over one period, -Rmax(^) = max{i?(t')|t < t' < t + T}. The origin 
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of the wiggles is a kind of resonance phenomenon in the RP equation and can quan- 
titatively be understood in detail Here we only discuss their consequences for 



the oscillating bubble. They mean that the bubble may stabilize through a tangent 
bifurcation: Imagine a fixed forcing Pa and then decrease Coo/cq. The Cqo/cq line will 
finally touch {p)^^/Po at a local maximum and create a pair of stable and unstable 
fixed points which will separate for decreasing Coo/cq. This process can repeat many 
times. The fixed points vanish through inverse tangent bifurcations, i.e., stable and 
unstable fixed points merge. The full bifurcation diagram is shown in figure ^ for 
Pa = l.latm and for Pa = l.Satm. 

The interpretation of the Rq - Coo phase diagram is as follows. The line signals 
equilibrium for the given driving pressure Pa- It is stable if the slope dR^/dcoolPa is 
positive, and unstable if it is negative. To the left of the line the bubbles shrink, to 
the right of it they grow by rectified diffusion. The shrinking or growing bubbles can 
hit a stable fixed point and thus stabilize. The basin of attraction of the stable fixed 
points is considerably larger for larger forcing Pa, see figure ||. 

A similar looking phase diagram results when fixing Cqo/cq and varying Pa- In 
figure 1^ we show such a bifurcation diagram in the Rq — Pa parameter space for large 
concentration Cqo/cq = 0.5, for lower concentration Coo/cq = 0.02, and for very low 
concentration Coo/cq = 0.002. Again, the lines signal equilibrium; stable equilibrium 
for positive slope dR^/dPalc^, unstable equihbrium for negative slope. To the right 
of the line we have growth, left of it shrinking. Increasing Pa at fixed c^o again leads 
to stabilization through a series of tangent bifurcations and later on destabilization 
by inverse tangent bifurcations. 

From these figures ^ and ^ we immediately understand why there is no diffusively 
stable SL for large At concentration [^]. No diffusively stable bubble radii exist in 
the large Pa - small Rq parameter regime where the energy focusing condition for SL 
is fulfilled. Note that there are stable equilibria, but for smaller Pa and large 
Rq- Note also that the basin of attraction of these stable equilibria is tiny. 

For small concentration Coo/ cq the situation is quite different. As seen from figure H 
and ^ now there are stable equilibria in the high Pa - low Rq regime where the bubble 
is both stable towards shape oscillations and fulfills the energy focusing criterion. 

In order to predict which SL regime is realized for given Cqo and Pa, we also have 
to take above shape instabilities into consideration. Therefore, in figure |10| we plot 
an enlargement of figure |^ together with the thresholds for the shape instabilities 
and the energy focusing condition, taken from figure |[ In the relevant parameter 
regime Pa > l.Oatm the critical R^^ is about 4 — 5/im. Only for large Pa > 1.4atm 
the Rayleigh- Taylor instability becomes relevant, see below for the consequences. As 
worked out in sections H and |T|, SL can only occur in the shaded region in figure |10[ 
In this region only for very low concentration Cqo/cq ~ 0.002 the bubble is diffusively 
stable and only then we can have diffusively stable SL. For larger concentration we 
have diffusively unstable SL. 

The different regimes of diffusively stable and unstable SL and no SL can easily 
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be read off from figure Talce fixed Pa = 1.1: growing bubbles are possible for 
Coo/ Co > 0.07, if they are large enough. They finally run into the parametric shape 
instability at ~ 5/im and a micro bubble will pinch off. If the remaining bubble is 
still large enough, i.e., above the stability line in figure H, the process will repeat. The 
allowed size of the bubble after the pinch off is very much restricted; for Coq/cq = 0.1 
it is only in the Rq range between 4.6/im and 5/im. If the pinched off micro bubble 
is larger, the remaining bubble will dissolve. No diffusively stable regime exists. 

For Pa = 1.3 the situation is quite different. For concentrations in a window 
0.003 < Coo/ Co < 0.005 bubbles in the Rq regime ~ 2fim up to 5/im will grow and 
stabilize. As in this Pa regime the energy focusing criterion is fulfilled, we have stable 
SL. Here we are at the very core of Barber's finding ||^ that stable SL is only possible 
for low Coo and high Pa- If the relative concentration is even smaller than 0.003, 
bubbles will dissolve and no SL is possible. For larger concentration large enough 
bubbles > l.Sfim (almost independent of Coo) will grow up to ~ 5/im where micro 
bubbles pinch off. In contrast to smaller Pa ~ l.latm, the remaining bubble is very 
likely to end up in the (now much larger) growing Rq regime. Pinched off micro 
bubbles with Rq < 1.8/im will dissolve. 

The same domains no SL, unstable SL, and stable SL can be identified from figure 
^ or figure 0; also these plots make it very evident that high Pa/ Po and low Coo/cq 
are necessary to obtain stable SL. 

The total phase diagram in the Coo - Pa phase space, our main result, was already 
presented in figure |l|. The notation in that diagram is as follows: If we denote a 
regime with stable SL or unstable SL we mean that there are bubbles of certain 
radius which are diffusively stable or growing, respectively; other, smaller bubbles 
dissolve. For Pa < 1.17 the Cqo/cq window of stable SL shrinks to zero. For very 
low Coo the no-bubble regime is very extended. If one now slightly increases Cqo, one 
immediately enters a regime where the energy focusing condition is fulfilled. This 
may be the prime reason why it is so much easier to find SL, diffusively stable or not, 
for low concentration. 

The experimental observations of fig. 4 in ref. |^ are in agreement with our anal- 
ysis. For that figure Pa is in the range of 1.3atm. Then we have unstable SL for 
large concentrations Cqo/cq = 26% and 6.6% and stable SL for low concentration 
Coo/Co = 0.4%. 

What happens for very large forcing amplitudes Pa? From figure 110 we see that the 
SL regime (shaded) intersects with the regime where bubbles cannot survive the micro 
bubble pinchoff (right of the dotted line (p. 18 ) in figure |10|) because of the violent 
short time scale Rayleigh- Taylor instability. Consequently, in this high Pa regime 
only stable SL should be possible, a prediction which is worth while being tested 
experimentally. Extremely low argon concentrations Cqo/cq < 0.0006 are necessary to 
guarantee diffusively stable equilibria in this regime. According to the approximate 
curve ( |2.18| ) of the onset of the RT instability in figure 0, this regime starts at about 
Pa ~ 1.4atm, but as pointed out in subsection II-D, this number should not be taken 
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too strictly. For lucidity we did not put in the right borderline of the unstable SL 
regime in the phase diagram figure |l|. - The upper borderline of the unstable SL 
regime is discussed in subsection IV-E. 

In our former publication we have stressed the importance of wiggles such as 
those in the graphs in figures |^ - As a consequence, for large enough Pa and small 
enough Rq one may have stable SL with different, discrete Rq; i.e., multiple stable 
equilibria. We will present some results of full numerical simulations on this issue in 
section 0. 

Possibly an experimental hint to multiple stable equilibria has been found by 
Crum and Cordry After registering SL light from a bubble for a few seconds, 
they distorted the bubble and observed the SL intensity to jump from one constant 
value to another (smaller) one. They suggested that the discrete light intensities 
corresponded to discrete diffusively stable radii. 

For the dynamics chosen here the multiple stable equilibria are just parametrically 
unstable as they lie beyond Rq ~ 5/im. This should not be taken too strictly, as 
a slight change of the model may allow for observable multiple stable equilibria. 
However, we see that the wiggles are not necessary for stabilization. An increase of 
(p)^4 (Rq) with Rq is sufficient, i.e., /3 > 0. If there are wiggles, we in addition have 
discretization of the equilibria. 



D. Scaling laws 

The average slope and the wiggles in figure ^ can be understood in detail as shown 



in ref. |^^. Here we only quote scaling laws for the running average of {p)^ ^ (Rq) which 



smooths out any wiggles. 



1 rRo+ARo 

^"^''-lAR^L.^ + (4.13) 



with ARq > O.Syum. For very small Rq the surface tension term dominates in the RP 
equation ( |1.2| ) and {p)^ ^ oc Rq^. For slightly larger _Ro < Rq"^ we have {p)^ ^ oc Rq^^'^ ■ 
For large Rq > Rq"^ (where wiggles occur and the running average becomes necessary) 
the first and the second term on the rhs of ( |L2| ) balance and 

W)t,A RT, (4.14) 

i.e., the average slope f3 is positive and equilibria are stable. The critical ambient 
radius Rq^^^ beyond which wiggles occur and stability is achieved scales like 



« 7^^, (4.15) 



the corresponding average pressure like 
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The scaling law ( [4.14|) directly reflects the average slope Rq oc c^^ in figure |[ 
Thus for large enough Rq there is always diffusive stabilization, but only for low Coo 
and high Pa this will be in a parametrically stable domain where the energy focusing 
condition is fulfilled. 

A detailed discussion of these types of scaling laws will be presented elsewhere 



E. Diffusively unstable SL bubbles 

In the previous subsections we saw that our theoretical stability diagram is in 
agreement with experiment. What about the growth rates in the unstable SL regime 
where the bubble grows by rectified diffusion and finally hits the parametric instability 
line at about Rq = 5fim7 Because of the growth (i) the relative phase of light emission 
0(i?o) will slightly change, (ii) the light intensity will increase as more and more gas is 
in the bubble, and (iii) the maximal radius will increase. When the shape instability 
line is hit, a micro bubble pinches off, giving the bubble a recoil. As this repeats 
again and again on the diffusive time scale of ~ Is, the bubble seems to "dance". 

Within the adiabatic approximation we now calculate ip^Roit)) = (pit). According 
to the energy focusing mechanisms discussed in section |T| light can be emitted if the 
(inward) bubble wall velocity becomes supersonic, M = —R/cgas > 1- We define the 
time when M = 1 holds as tg- The relative shift of this time to the forcing phase 
defines 0s . As the waves and shock waves in the bubble are very fast, we take the 
time tg of the detachment of the (shock) wave as the time of the light pulse. The 
error we make by this approximation is of the order of At ~ Ro/cgas ~ Ins, as follows 
from a simple estimation. 

Next we calculate -Ro(^ from eq. ( [4.1(]| ) for discrete adiabatic times t = nT. Take 
fixed Rq and calculate the time averages {.)^- from the RP equation. The integral 

(4.17) 



m+R%t)m,, 

can be calculated numerically. Its convergence, however, is slow. In order to speed up 
the numerical calculation, for calculations over a long period of time we approximate 
the integral / by / = a/Rmax + (1 — 0')/Ro + (3/imax)^^^'^- Here, a ~ 0.9 is an 
adjustable parameter, which slightly depends on Rq and Pa, and hmax ^ Rmax must 
be sufficiently large. The approximation is very well controlled and the results are 
indistinguishable from the exact result. The growth during the time interval T finally 
reads 



Co Pq 



(4.18) 
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ARq is added to Rq and the procedure is repeated until -Ro(^ hits the parametric 
instabiUty curve. Here from figure ^ we took R^^ = 5/im as a very good approxi- 
mation. A random fraction of the bubble will pinch off. In figure |Tl| we show the 
ambient radius -Ro(^ and the relative phase of the light pulse 4>s(t) for three different 
relative gas concentrations Coo/cq. 

Figure [l^ should be compared to the corresponding experimental figure 4 of ref. 

0. Unfortunately, for that figure the precise forcing pressure amplitude Pa and the 
ambient bubble size are not known. So we try both Pa = 1.3atm (fig. |ll|a) and Pa = 
1.2atm (fig. lib). For Pa = l.Satm we have very good agreement with experiment: 
stable SL for 3mmHg, growth for SOmmHg and for 200mmHg. For Pa = 1.2atm a 
bubble in a fluid with 3mmHg argon concentration would dissolve according to our 
approximation, see figure |I], so we choose 12mmHg as smallest concentration and find 
stable SL. 

Let us compare the growth rates with the experimental values. As in experiment, 
for larger gas concentration the growth rates of Rq and (pg strongly increase. Quan- 
titative estimates from figure |n| for the growth rates of the phases for both forcing 
pressures and both concentrations are compared to the experimental ones in table 

1. They are slightly larger from what is found in experiment, but agree in order 
of magnitude. Rather than focusing on an exact quantitative agreement here, the 
important point is that for lower concentration Cqo/cq = 0.0165 (corresponding to 
3mmHg) the phase of light emission is stable due to the diffusively stable ambient 
radius (for Pa = l.Satm). In this example, the bubble is locked at Rq = 4.31/im, and 
the phase of the light pulse is stable for ever. 

What is the physical consequence of the large bubble growth rates obtained for 
large argon concentration Cqo/cq (table 1)? The shape instability line will be hit more 
frequently per unit time and the bubble's dancing frequency will thus become larger, 
as more micro bubble pinch-offs and resulting bubble recoils will take place per unit 
time. As noted above with some certain probability per pinch-off the pinched off 
micro bubble (s) are too large so that the remaining bubble dissolves. Thus with 
increased pinch off frequency this probability per time increases. We speculate that 
this mechanism sets the upper threshold of the unstable SL regime (towards a no 
SL regime) in the phase diagram figure |I[ Table 1 teaches us that the growth rates 
increase drastically with Coo/cq, so the probability of having a long living bubble for 
high gas concentration becomes very low. Indeed, Gompf reported that the larger 
the concentration is for fixed Pa, the faster the unstable SL bubble dies. If Coo is big 
enough, it will thus be very unlikely for the bubble to survive an appreciable time. 
From experiment we know that the upper concentration threshold of unstable SL 
is beyond Cqo/cq = 0.26. 

This hypothesis also explains why the water in the SL container "ages" in 
case the container is not gas tight. By "aging" it is meant that stable SL and finally 
also unstable SL becomes impossible with "old" water. The reason is that external air 
diffuses into the water, dissolves, and Cqo/cq increases. Consequently, the originally 
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stable bubble is pushed into the unstable regime and starts to "dance", shedding off 
micro bubbles. The dancing frequency becomes higher and higher and finally the 
bubble dissolves after a too large pinch off. Bubbles may be reseeded, but will also 
die very soon. 

What we do not understand in the unstable regime is the dependence of the light 
intensity on the gas concentration as e.g. measured in figure 2 of |T^. We speculate 



that it depends on the ambient size of the bubble which is supported by figure 6 of ref. 
0. In that figure Lofstedt et al. show that the maximal radius and the SL intensity 
are correlated. 





Coo 


Pa = 1.2atm 


Pa = l.Satm 


experiment 




SOmmHg 
200mmHg 


1.7 fis/s 
9/is/s 


2.7 fis/s 


0.5/is/s 

5fis/s 


Ro 


SOmmHg 
200mmHg 


2.4/im/s 
llfim/s 


4.7 fim/ s 
18/im/s 





Table 1 

Growth rates for the phase of light emission (ps for two forcing pressure amplitudes Pa 
and two concentrations Coo near the pinch off of the micro bubble in comparison with 
the experimental data /§/ for which the forcing is not exactly known. Stronger forcing 
and larger argon concentration enhance the growth. We find order of magnitude 
agreement. We also give the growth rates for the radii which are not experimentally 
available. 



Finally we address the question whether the growth of 0s has to be monotonous. 
In fact, it does not. For larger Rq Gfim, where the bubble is already parametrically 
unstable) the growth rate of (pg is wiggly. Similar oscillations show up in the maximal 
radius -Rmax(^) = max{i?(t')|t < t' < t + T}, as e.g. seen in figure 4 of Barber et 
al. 10. This wiggly structure as a function of time is a direct consequence of the 
wiggly structure of 0^ and Rmax as a function of ambient radius Rq which is due to a 
resonance phenomenon in the RP equation The growing Rq probes the wiggles 
in (j)s{Ro) and Rmax{Ro)- 

In experiment no or hardly any oscillatory structure in 0^ is seen 0. Consequently, 
the ambient radius Rq does not seem to be in the wiggly regime. Thus multiple stable 
diffusive equilibria may only be important in the shape unstable regime in the Ro — Pa 
parameter space. Indeed, the wiggly structure in Rmax{Ro(t)) is only revealed when 
boosting the bubble in the unstable SL regime where it becomes shape unstable after 
a few ms and bursts (cf. fig. 4 of 0). 
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V. COMPARISON OF THE ADIABATIC APPROXIMATION TO THE 

FULL NUMERICAL SOLUTION 



A. Concentration profiles 

To compare our results within the adiabatic approximation with the exact solu- 
tion, we must numerically solve the PDE ( |4.1| ) with the boundary conditions given 
above. We sketch our numerical method in appendix ^ Here we report on results. 
We do not consider shape instabilities in this section. 

In fig. we show concentration profiles of the gas outside the bubble during 
expansion and collapse. We take an argon bubble driven at l.lSatm close to diffusive 
equilibrium. Near the minimum the concentration gradient at the bubble wall is 
negative and the bubble ejects gas which accumulates near the bubble wall as the 
diffusive time scale is slow compared to the bubble motion. When the bubble is 
reexpanding, it pushes away the accumulated gas together with the fluid. During 
the expansion phase the concentration c{R{t),t) at the bubble wall decreases due to 
Henry's law ( [4.2|) . At some point the gradient becomes positive and the mass content 
of the bubble grows. The wall of gas outside the bubble is thus (i) pushed away from 
the bubble, (ii) deaccumulates (towards smaller r) because of bubble growth, and (iii) 
shrinks because of diffusion (towards larger r). 

Apart from the diffusive processes in the fluid these proflles only mirror the time 
and space dependence of Cosc{r,t). The width of the boundary layer in which Cosc is 
a dominant feature of the concentration proflle can be readily estimated as 

Sd = (5.1) 
V 00 

which is the characteristic length for diffusive processes on a time scale ~ T. 



B. Comparison on diffusive equilibria 

First, let us focus on the equilibrium radii for diffusively stable bubbles. Can we 
flnd the multiple equilibria predicted in section |^ in the full numerical simulation, i.e. 



different stable situations for the same physical parameters? Figure |13| shows Ro{t) 
for several different initial ambient radii. The diffusion constant has the physical 
value D = 2 ■ 10~^m^/s, corresponding to a Schmidt number of Sc ~ 3000 in the 
regime of interest (cj ^ 27r • 26.5kHz, Rq ~ Gfim). We indeed observe several stable 
and unstable equilibria. However, they deviate slightly from those calculated in the 
adiabatic approximation. 



The deviation, however, is tiny, less than O.OG^m as seen from flgure and clearly 
not detectable with today's experimental possibilities. Thus for all practical reasons 
we can consider the adiabatic approximation of the equilibrium radii as exact. 
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Nevertheless, let us wonder what the origin of the deviation is. It can be explained 
by considering higher order corrections to the adiabatic theory [^: the equilibrium 
condition ( [4.11| ) is modified, lowering the required mean pressure at the bubble surface 



for stability to 



W))tA Coo 1 4/1 



Po Co 5ci/2 Co \R^Jo ' ' I 



c. 







.{h,t)dh) , (5.2) 



tA 



Cq^c being the zeroth order (in Sc''^^'^) solution of the oscillatory part of the profile, 
depending on h = Sc^^'^h IQ. The shift in {p{t))^^ causes a corresponding shift 
j^adiab _ j^e ^^le equilibrium radius R^. 

To analyze the deviations from the equilibrium position (|4.11| ) further, we redid 
the numerical PDE calculation for smaller (unphysical) D. As expected, in the D — > 
limit the {D independent) adiabatic fixed point is approached as can be seen from 



fig. |TJ. For that figure the numerical equilibrium radii are "measured" by detecting 
either "shrinking" or "growth" for slightly different radii. The stable equilibrium 
is determined by linear interpolation between the growth rates of a growing and a 
shrinking bubble. In principle, we can achieve arbitrarily small error bars in our 
numerical results. This is confirmed by the excellent agreement of the D — > limit 
with the adiabatic approximation. 

There is a theoretical possibility that the deviation of the full PDE dynamics from 
the adiabatic approximation matters, namely when multiple stable equilibria are to 
be resolved. For completeness we discuss this point in Appendix O. 



C. Comparison on growth rates and profiles 

To detect the location of an equilibrium point along the Rq axis, starting with a 
constant gas concentration profile c{h) = Coo is obviously a good choice, because Cgmo 
will approach that value for all h in equilibrium (apart from higher order corrections). 
The sign of Ro{t) will be correct after a small number of oscillation cycles. If one is, 
however, interested in the actual value of the bubble growth rate, i.e. bubble dynamics 
far from equilibrium points, choosing Csmo{h) as initial concentration profile will avoid 
transients on diffusive time scales. Indeed, 

{c{h,t))^^ ^ Csmo{h) (5.3) 



holds to very good accuracy for such a calculation (see figure ITS]). Moreover, the 
observed growth rate -Ro(^) is in very good agreement with the value calculated from 
( [4.101 ) , as indicated in figure |T6|. 
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VI. CONCLUSIONS 



We worked out a hydrodynamic approach towards SL, based on the Rayleigh-Plesset 
equation. This allows us to explore a considerable part of the phase space and to study 
long term dynamics. As necessary conditions for SL we demanded shape stability and 
energy focusing. The adiabatic approximation p7|,|7| allows us to study also diffusive 
stability within the RP approach. 

We have presented phase diagrams in the Coo — Pa, Ro — Pa, and Rq — Coo parameter 
spaces. Three phases can be identified: stable SL, unstable SL and no SL. Stable SL 
only occurs in a tiny domain of the parameter space which is in good quantitative 
agreement with the UCLA SL experiments on argon bubbles |^. 

For molecular gases besides (i) shape instabilities and (ii) diffusive instability also 
(iii) chemical instabilities have to be considered. Then our hydrodynamic approach 
can be extended to gas mixtures such as air as demonstrated in refs. [^,^. Again, 
good agreement with the UCLA experiments [|,^^|^,|TT| is achieved. 

We suggest to experimentally map out phase diagrams in the Cqo/cq versus Pa/ Pq 
parameter space for various gas mixtures. Detailed experimental data on the border- 
lines between the stable, unstable, and no SL regimes will lead to further refinements 
and improvements in our understanding of the dynamics of a sonoluminescing bubble. 

The central question which cannot be answered within the present approach is how 
hot the gas inside the bubble can become. Progress on this point will require more so- 
phisticated understanding of the gas dynamics inside the bubble. Understanding how 
the gas temperature depends on experimental parameters such as forcing pressure, 
gas concentration, or liquid temperature will allow for the creation of temperature 
controlled environments for chemical reactions within the bubble. 
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APPENDIX A: NUMERICAL SCHEME FOR THE ADVECTION 

DIFFUSION EQUATION 



1. Transformation 



To numerically solve eq. (|4.1| ) with the boundary conditions ( [4 .21) and (|4.3| ) we 
first transform the independent space variable r. Following we choose 



F{x,t) = {{R{t))'-Rl\n{l- 



1/3 



(A.l) 



where Ru is an adjustable length parameter. All lengths are measured as multiples 
of lyum, all times in units of T = 27i/uj, and all pressures in atm. After the transfor- 
mation r = F{x, t) the advection diffusion equation ( [4. ID reads, if written as a local 
conservation law, 



dtifc) + d.^{gc) - d^{hd^c) = 0, 



(A.2) 



where 



f{x,t) 
g{x,t) 
h{x, t) 



F^d^F. 
RR — ±- ut 
DF^/d^F. 



2 - F^aF, 



(A.3) 



The transformation ( |A.1| ) is chosen to make g{x,t) identically vanish and thus to 
obtain a pure diffusion equation. We have 



fix,t) 
g{x,t) 
h{x, t) 



0, 

D{l-x){F{x,t)r/Rl. 



(A.4) 



The time dependent range r G [R{t), oo] has been mapped to the constant range 
X G [0,1]. The boundary conditions are Henry's law ( [4.2| ) which now reads c{x = 
0,t) = cop{R,t)/ Pq and c(x = l,t) = Coo- The mass loss of the bubble, expressed in 
the new independent variable x, is 



m 



AnDR-^R'^drC 



(A.5) 



2. Discretization 

The interval x G [0, 1] is discretized using a non-equidistant grid. The grid must 
(i) provide sufficient volume for diffusion of gas outside the bubble and (ii) sufficient 
resolution near the bubble radius for a correct representation of gas concentration 
gradients. We satisfy (i) by choosing _R„ ^ Rmax, where Rmax is the maximum radius 
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of a typical bubble oscillation. This ensures r^{x ^ 1) ^ even for relatively low 



resolution near x = 1 (cf. eq. (|A.1|) ). From this choice, we deduce a criterion for 
the grid resolution Sxr near the bubble surface for a given resolution Srn in physical 
coordinates (e.g. Srn = O.Ol/xm). Expanding ( |A.1|) in x -C 1, we get 



6x 



R 



max 



This is a conservative estimate, requiring 5r/j to be resolved even for the minimum 
bubble radius Rmin- At times when the bubble radius is larger, the resolution in r 
will be even better with this Sxr. 

Typical bubble dynamics data lead to 6xji ~ 5 • 10^^. Excellent resolution at the 
bubble radius is needed for a correct representation of gas concentration gradients. It 
is, however, unnecessary at greater distance from the bubble surface. Therefore, the 
grid density is varied according to a power law to yield a fine grid near a; = and a 
coarse grid near x = 1. 

Overall, we have grid points Xi = 0, X2, . . ., xjy^i, xn = 1. The field c{x,t) is 
defined on the grid, Cj = c(xj), and so are the fields /, g, h. We define dxi = Xj+i — Xi, 
i = 1,2, ... N — 1 and dxa,i = {dxi-i + dxi)/2, i = 2, . . . N — 1, dxa^i = dxi/2, 

dXa,N = dXN-l/2. 



The discretization of eq. (|A.2|) has to guarantee mass conservation. 











c{-piR%t))co/Po, 

•1^ dt dxai \ 2 V^i ^ "-i+lf dx. 



2....N -2. 



> (A.6) 



dcj 
dt 



AT-l, 



Here, are the concentrations q + Odci at time t + 9dt and /if = hi{R^) where R^ is 
the radius at time t + 9dt. Correspondingly, Rq is the ambient radius at that time. We 
choose = 1, i.e., a fully implicit method. Equation ( |A.6|) has to be assisted by the 
Rayleigh-Plesset equation ( |1.2| ) and the proper discretization of ( |A.5| ) guaranteeing 
total mass conservation, 







hi + hi 4 



dxi 



(A.7) 



We solve the A^ + 2 equations ( |A.6D , ( |A.7D , and (|1.2|) for the unknowns dci, i = 
1,2, ... N, dRo, and dR with Newton's method. The Jacobian is calculated analyti- 
cally. A time step control and adjustment is provided by redoing every time step dt 
in two steps of width dt/2 each and then comparing the result. For large forcing we 
need a very low tolerance of 10~^ per cent to achieve sufficient numerical quality of 
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mass conservation. Note that this simulation covers time scales from picoseconds (for 
good resolution of the bubble dynamics near the collapse) to tens of milliseconds (for 
observation of diffusive growth or shrinking of the bubble). 



APPENDIX B: ADIABATIC APPROXIMATION AND MULTIPLE 

EQUILIBRIA 

In section we showed that the adiabatic approximation well describes the full 
dynamics. However, for (theoretical) completeness we would like to caution in this 
appendix: For the physical diffusion coefficient D = 2- 10~^m?/ s the simple adiabatic 
approximation may lead to even qualitatively wrong results. E.g., take a bubble with 
Ro = 6.06/im. According to the adiabatic approximation one would expect that it 
grows towards the equilibrium R'^'^^°-^ = 6.097/im (see figure 0). However, it shrinks 
towards the physical equilibrium Rq = 6.041/im. At the moment the experimental 
accuracy does not allow to distinguish between these two sizes, but it may improve 
one day. We can understand Rq < R^'^^"-^ because the gas layer around the bubble 
(figure p!^ ) diffusively shrinks too strongly for finite D, leading to a larger overall 
mass loss and thus smaller bubbles. The order of magnitude of the R^ shift can be 
estimated as the width of the boundary layer: 

jadiab rye 




Rf'^''-Rl^6D = \-ocD^/\ (B.l) 



For the physical D = 2 ■ lO^^m^/s we have R'^'^'^"-^ — Rq ^ O.ll/im reproducing the 
numerical result ~ 0.06/im quite accurately. In the inset of fig. |l^we plot log{RQ'^'^°'^ — 
Rq) vs logD and indeed find good agreement with the scaling law ( [B . 1| ) . 

Whether the deviations are considered to be serious or not depends on what is 
supposed to be analyzed. If the focus lies on identifying equilibrium points, the con- 
sequences can be quite drastic. Note that for large enough D, the value of (|5^) falls 
below the minimum values of {p{t))^^/PQ (cf. figure]^), thus making an equilibrium 
solution impossible and leading to dissolution of the bubble. 

It the distance between subsequent equihbrmm radn Aito =-".0 Rq 
(where the index numbers label the equilibria in ascending Rq order) is supposed to be 
resolved, i.e., if one wants to find a one to one correspondence between the adiabatic 
equilibria and the real ones, one has as to impose the following condition on the 
approximation: 

^j^adiab,i ^ ^^j^adiab,i _ ^e) _ 2) 

Expressed in the Schmidt number Sc it reads 

(B.3) 
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rather than simply 5*0 S> 1 as one may naively expect. With the correct values in the 
relevant parameter regime Ro ~ 6/im, ~ O.Sfim we have Sc ^ 1600, which 

is only marginally fulfilled by the physical D corresponding to Sc a; 3000. 
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FIG. 3. Time development of the normalized distortion amplitude a2{t)/R{t) for (a) a 
Rayleigli- Taylor unstable, parametrically stable bubble (i?o = 2.5/2771, Pa = l.Satm), (b) an 
afterbounce unstable, parametrically stable bubble (i?o = 4.0/Ltm, Pa = l.Satm) and (c) a 
Raylcigh- Taylor stable, parametrically unstable bubble {Rq = 5.2//m, Pa = l.Oatm). In (a) 
wc also show a2/R in a blow up of the time scale to demonstrate that the typical time scale 
of the Rayleigh- Taylor instability is nanoseconds. The typical time scales of afterbounce (b) 
and parametric instability (c) is in the microsecond, and millisecond range, respectively. 
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FIG. 4. Time development of the bubble radius R{t) (upper part) and distortion ampli- 
tude 02(t) (lower part) for a i?o = ^Ajim bubble driven at Pa = l.latm. Note the transition 

from Rayleigh- Taylor (time scale ns) to afterbounce perturbations (time scale jis) during 
the afterbouncc part of the bubble dynamics. It is also seen that the dynamics of the dis- 
tortion a2{t) has half the frequency of its forcing bubble dynamics R{t) as typical for an 
instability of the Mathieu type. 
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FIG. 5. Borderline of the parametric instability (solid), the afterbounce instability 
according to criterion ( |2.15| ) (short dashed), and the M = 1 criterion (|3.1| ) for a supersonic 
bubble collapse (long dashed). The Re = 10- criterion ( |3.3| ) for the persistence of a shock 
(dot-dashed) is found to be less stringent than the M = 1 criterion. Also shown is the 
perturbation velocity threshold (2.18) (dotted). In case of unstable SL, to the right of this 
curve the bubble is thrown out of the trap as the pinch off of micro bubbles is too violent. 
The region where SL is possible is shaded. 
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FIG. 6. Maximum perturbation velocity according to ( p.l8D as a function of Pa and Rq. 
This calculation was done with an 02 noise amplitude of Inm. Note the two distinct regions 
of the parametric instability (left) with its rapid onset and the Rayleigh- Taylor instability 
(right, high pressures and small Rq). In the latter region the bubble is thrown out of the 
trap when the pinch off occurs. It corresponds to the region right of the dotted line in figure 
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FIG. 7. {p)t,A/PQ as a function of the ambient radius Rq for small forcing pressure 
amplitudes Pa=0-8atm to Pa=l-15atm (left, top to bottom, in steps of 0.05 atm) and for 
large Pa=1.2, 1.25, 1.3, 1.4, 1.5atm (right, top to bottom). Note the different ordinate 
scales. 
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FIG. 8. Bifurcation diagrams in the Rq — Cqo parameter plane for a forcing pressure of 
Pa = 1.3atm (left) and Pa = l.latm (right). Tangent bifurcations are seen. The regimes 
with positive slope are stable. To the left of the curves the bubbles shrink, to the right of 
them they grow. 
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FIG. 9. Bifurcation diagrams in the Rq — Pa parameter space. Tangent bifurcations are 
seen. The regimes with positive slope are stable. Gas concentrations are Cqo/cq = 0.002 
(right), Coo/co = 0.02 (middle), and Cqo/cq = 0.5 (left). To the left of the curves the bubbles 
shrink and finally dissolve, to the right of them they grow by rectified diffusion. 
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FIG. 10. The figure shows the discussed effects for argon all together: Beyond the M = 1 
curve (long dashed) SL is possible. The bubble grows thanks to rectified diffusion right of 
the diffusive stability curves (shown for Cqo/cq = 0.5, 0.02, and 0.002, left to right). At the 
parametric shape instability (solid) and the shape instability according to criterion ( 2.15 ) 
(dashed) micro bubbles pinch off, but the bubble can survive in the trap for Pa < 1.4atm. 
SL is possible in the shaded region, stable SL, if in addition the slope of the diffusive 
equilibrium curve is positive. Here, this occurs for Coo/co = 0.002 for Pa ~ 1.33atm and 
-Ro ~ 3//m. Right of the dotted line (criterion (|2.18|) ) the bubble is thrown out of the trap in 
case of a pinchoff. Consequently, if at all, only stable SL is possible in this high Pa regime. 
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FIG. 12. Gas concentration profiles outside the bubble for expanding (a) and collapsing 
(b) bubble radius, respectively. The bubble is near diffusive equilibrium, driving pressure 
is Pa = l.lSatm. Profiles are shown at intervals of 0.75/xs (a) and 10ns (b). 
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FIG. 13. Three different stable equilibria are approached both from above and from 
below. Shrinking or growth are hardly noticable. The six initial radii are 6.035;um, 6.05//m; 
6.32//m, 6.34/im; 6.65/im and 6.665/im, respectively. Otherwise, all conditions are the same. 
Again, we choose Pa = l.lSatm. 
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FIG. 14. Location of an equilibrium radius for various D (solid line). At the upper end 
of the error bar the bubble is shrinking, at the lower end it is growing. The theoretical value 
from adiabatic approximation is indicated by the dashed line. Note that the tiny difference 
of ~ 60nm (for the physical diffusion coefficient D = 2 ■ 10~^m^/s) is not experimentally 
detectable at the moment. The bubble is driven at l.lSatm with a gas saturation in water 
of Coo/co = 0.035. The inset shows the deviation R^'^^"-^ — drawn in a log-log diagram, 
together with our estimate (|B.1| ) (solid line). 
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FIG. 15. Comparison of (c(r, i))t,4 (solid line) and the smooth profile CgmoiT-, t) from the 
adiabatic approximation (dashed). The bubble ambient radius shows pronounced growth 
in this case. The inset shows an enlargement of the small r regime. 
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FIG. 16. -Ro(^) for Coo/cq = 0.001 and for Coo/cq = 0.3 resulting in shrinking and growth 
of the bubble respectively. We chose Pa = 1.15atm and a bubble with Ro{t = 0) = 5.5/im. 
The adiabatic approximations are also shown as straight lines. By definition they do not 
follow the violent mass exchange processes during one cycle of time T. 
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